Quark mass dependence of the pseudoscalar hairpin vertex 



W. A. Bardeen a , E. Eichten a , and H.B. Thacker b 
a Fermi National Accelerator Laboratory, P.O. Box 500, Batavia, IL 60510 
b Department of Physics, University of Virginia, Charlottesville, VA 22901 

(Dated: May 20, 2004) 

In a recent investigation of chiral behavior in quenched lattice QCD, the flavor-singlet pseudoscalar 
"hairpin" vertex associated with the rf mass was studied for pion masses ranging from « 275 to 675 
MeV. Throughout this mass range, the quark-disconnected pseudoscalar correlator is well-described 
by a pure double-pion-pole diagram with a p 2 -independent mass insertion. The residue of the double 
pole was found to exhibit a significant quark mass dependence, evidenced by a negative slope of 
the effective mass insertion (m^) 2 as a function of m\. It has been observed by Sharpe that, 
with a consistent NLO calculation in quenched chiral perturbation theory, this mass dependence is 
uniquely predicted in terms of the single-pole coefficient a<s> and the Leutwyler parameter L5. Since 
aj. is found to be ss 0, the chiral slope of the double-pole residue determines a value for L 5 . This 
provides a consistency check between the chiral slope of the hairpin mass insertion and that of the 
pion decay constant. We investigate the consistency of these mass dependences in our Monte Carlo 
results at two values of lattice spacing. Within statistics, the slopes are found to be consistent with 
the Q\PT prediction, confirming that the observed negative slope of mf^ arises as an effect of the 
I/5 Leutwyler term. 

PACS numbers: 11.15.Ha, 11.30.Rd 



The predictions of chiral perturbation theory provide 
extensive tests of numerical calculations for light-quark 
amplitudes in lattice QCD. In a chiral Lagrangian frame- 
work, the quenched approximation introduces a variety 
of anomalous effects associated with the flavor singlet 
pseudoscalar rf meson. These quenched chiral loop ef- 
fects can serve as a useful diagnostic for lattice QCD 
calculations with light quarks. Successfully reproducing 
the quark mass dependences predicted by chiral pertur- 
bation theory provides an important test of the practical 
validity of any chiral lattice fermion method. Although 
the theoretical underpinnings of quenched chiral pertur- 
bation theory (QxPT) 0, 0] are less rigorous than for 
full QCD, recent lattice calculations 0,0, have exhib- 
ited good agreement with QxPT predictions in a number 
of amplitudes involving scalar and pseudoscalar meson 
correlators. Within the mass range studied (m, « 275 
to 675 MeV), these results provide strong support for 
the assertion that the effect of the chiral U(l) anomaly 
is well-represented by a p 2 -independent rf mass inser- 
tion, which can be treated perturbatively according to 
the rules of QxPT- 

A direct study of the time-dependence of the singlet 
pseudoscalar correlator shows that the hairpin insertion 
vertex is a pure rf mass term, oc m 2 ,^' 2 , with higher 
derivative terms being suppressed. The a$ parameter 
of a possible a$(<9 M 7/) 2 term is found to be consistent 
with and definitely << 1. Specifically, the correlator 
fits well at all time separations to a pure double-pion pole 
formula, with little evidence for either a single pole term 
(parametrized by ev$) or for any excited-state contamina- 
tion. An interesting feature of the results presented in Q 
is that the effective rf mass insertion (m^) 2 obtained 
from the residue of the double pole decreases significantly 
with increasing quark mass. The possibility of explaining 



this mass dependence as an effect of the p 2 -dependent cv$ 
term is ruled out by the shape of the time dependence of 
the hairpin correlator, which gives a$ ~ 0. 

A consistent explanation of the mass dependence of 
the double-pole residue has recently been proposed by 
Sharpe 6], who observed that the value of (m^/^) 2 re- 
ceives an order m 2 contribution from the L5 Leutwyler 
term. Sharpe also showed that any other possible terms 
of order m q or p 2 that could contribute to the mass de- 
pendence of the hairpin double pole residue can be re- 
moved by a field redefinition along with a redefinition 
of cv$. As a result, such terms do not introduce any 
new physical parameters, and therefore the lowest order 
quark mass dependence of the hairpin insertion should 
be given entirely by the cv$ and L5 terms. Furthermore, 
estimates of both a$ and L5 can be obtained from the 
hairpin analysis alone by first determining cv$ from the 
residue of the single-pole term in the hairpin correlator 
(which gives cv$ « 0), and then attributing the remain- 
ing mass dependence of the double-pole residue to the 
-L5 term. Since the L5 term is also responsible for the 
quark mass dependence of the axial-vector decay con- 
stant /a, we have two independent ways of estimating 
-L5 and hence a consistency check on the QxPT analy- 
sis of pseudoscalar correlators. Note also that, for equal 
quark masses, quenched chiral logs are not present in ei- 
ther /a or rriQ . Thus linear quark mass dependence is 
expected and is in fact observed in our data 0. Taking 
a$ = 0, QxPT predicts ja| that the slope of is equal 
in magnitude and opposite in sign from that of /a- Thus 
the product x j 'a should be independent of quark 

mass, and the negative slope of the mass insertion may 
be inferred from the well-known positive slope of j a- In 
this brief note, we examine the quantitative consistency 
of our data with this prediction. 
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TABLE I: Values for m^,mQ^, and /a for f3 = 5.7. Values are 
in GeV, using the rho mass scale a -1 = 1.12 GeV. Values of Ja 
include a perturbative renormalization factor of Za = 0.845. 
Values of tUq should be multiplied by a flavor factor of \/3 
for comparison with the physical rf mass. 





m n 


eff 


/a 


.1410 


.566(2) 


.314(11) 


.173(2) 


.1415 


.504(3) 


.329(11) 


.167(2) 


.1420 


.432(3) 


.345(11) 


.160(2) 


.1423 


.383(4) 


.354(11) 


.156(3) 


.1425 


.344(4) 


.360(11) 


.154(3) 


.1427 


.299(5) 


.361(12) 


.152(4) 


.1428 


.274(6) 


.361(12) 


.152(5) 



TABLE II: Values for m^m^ 7 , and f A for /3 = 5.9. Here 
a' 1 = 1.64 GeV, and Z A = 0.865. 



K 




m o 


fA 


.1382 


.674(5) 


.318(8) 


.172(1) 


.1385 


.620(5) 


.325(8) 


.166(1) 


.1388 


.563(7) 


.333(8) 


.160(2) 


.1391 


.499(7) 


.343(8) 


.154(2) 


.1394 


.428(8) 


.356(8) 


.147(2) 


.1397 


.334(8) 


.371(10) 


.139(3) 



First we summarize the relevant Monte Carlo calcu- 
lations from Ref. [3] and review the QxPT predic- 
tions for the mass dependence of mj,^ and Ja @- The 
quark-connected (valence) and quark-disconnected (hair- 
pin) pseudoscalar correlators are defined in momentum 
space by 
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FIG. 1: The ratios Rf (open squares), Rm (x's), and the 
product RfRm (solid squares) for (3 = 5.7 (upper plot) and 
P = 5.9 (lower plot). 



Ac(p) 



d 4 xe lp - x 



{qAi^qsix) gBj5qA(0)) (l) 



and 



d 4 xe tpx 



(lAl5qA{x) (?S75(?S(0)) (2) 



Here qA and qs are two quark fields with distinct fla- 
vor indices. For the results considered here, we take 
the two quarks to have degenerate mass. Clover im- 
proved Wilson-Dirac fermions are used, with C sw = 1.57 
for the (3 = 5.7 configurations and C sw = 1.50 for 
(3 = 5.9. Quark propagators are calculated in the mod- 
ified quenched approximation (MQA) by locating and 
shifting real eigenmodes, as described in |7j. Fitting the 
long-range time dependence of the Fourier transformed 
p = valence correlator in the usual way, we extract 
the pion mass and pseudoscalar decay constant from the 
location and residue of the propagator pole, 



Ac(p) 



Pp 



(3) 



The values obtained for are listed in Tables I and II, 
along with the values for the axial- vector decay constant 
Ja- The complete analysis is discussed in 0- We do not 
list the values of fp since they will not be used in the 
present analysis, except to amputate the corresponding 
factors in the hairpin correlator. With pion masses fixed 
from the valence correlator analysis, the hairpin correla- 
tors are fit to a combination of double- and single-pole 
terms, 



A h (p) = {ml ff f- 



(p 2 



-,2\2 



fp 



(4) 



The values obtained for m JJ are listed in Tables I and 
II. From the Q\PT analysis in Ref. 6], it is found that 
the double-pole residue is given, to first order in m^, in 
terms of the Lagrangian parameters mo, L5, and /, 
by 



/ e 

(to 



e/j\2 = 



1 8T m - 

1-8L 5 — 



(5) 
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As discussed in Q the single pole coefficient a$ is small 
and consistent with zero. In the subsequent analysis, we 
will assume a$ = 0. With this simplification, Q\PT pre- 
dicts that the chiral behavior of mjj^ is related to that of 
/a by the requirement that their product is independent 
of quark mass to first order, i.e. 



MS 



x Ja = const. + 0(m*) 



(6) 



To directly test this prediction, we would like not only to 
see that the slope of the product is consistent with zero, 
but also that a significant slope in each of the two quan- 
tities separately is largely cancelled out in the product. 
This is conveniently illustrated by plotting the ratio of 
each quantity to it's chirally extrapolated value. Define 
the quantities 



RsiO 

R m (ml) 



UK) 

/a(0) ' 

w 



in 



efS 



(7) 
(8) 



In Figure 1 we have plotted the quantities Rf, R m , and 
the product Rf x R m as a function of m \ . The errors that 
are plotted were computed by a jackknife analysis, with 
chiral extrapolation carried out on each subensemble to 
determine the denominators in (7J| and (JSJ). The plots 
in Fig. 1 make it quite clear that the slopes of Ja and 
m|j^ approximately cancel in the product. For the data 
at P = 5.7 the cancellation is quite accurate and the 
slope of the product is very close to zero. For (3 = 5.9, 
the product appears to have a slight positive slope. A 
jackknife analysis directly on the slope of the product 
gives 



Rf Rn 



1 



em, 



(9) 



with e = -0.051 ± 0.16 for (3 = 5.7 and e = 0.20 ± 0.14 
for (3 = 5.9 to be compared with the individual slopes of 
about ±0.60 for Rf and R m separately. The difference of 
this slope from zero in the /3 = 5.9 ensemble is of marginal 
significance and may be purely a statistical fluctuation. 
It should be noted, however, that the numerical value of 
the clover coefficient C sw for (3 < 6.0 is not very well 
determined by theoretical considerations. It may be that 
our choice of C sw = 1.50 for (3 = 5.9 leaves some residual 
0(a) lattice spacing effects. The sensitivity of the slope 
in (0 to choice of clover coefficient has not been explored. 

As emphasized by Sharpe [||, the quark- mass depen- 
dence of the hairpin vertex provides a new way of de- 
termining the quenched Leutwyler parameter L§. The 
extraction of the Leutwyler parameter L§ from the slope 
of Ja was discussed in j3|. By a similar analysis, an es- 
timate of L5 may be obtained from the slope of m e . 



Although the difference between these two determina- 
tions is not very significant statistically, the values ob- 
tained from the m[j^ slope are somewhat smaller than 
those obtained from /a, particularly for the 5.9 data. 
From the jackknife subensemble analysis, we find that 
the slopes of Ja and are positively correlated, i.e., 

since the slopes are of opposite sign, the magnitude of 
the two slopes (and hence the corresponding estimates of 
i5)fluctuate in opposite directions. This suggests that we 
compare the slope of Ja directly with the slope if l/rn^^ . 

To this end, we analyze the slope of Ja, \J j ' Afm^ ^ and 

I/tUq in the form A + Brri^. The Leutwyler parameter 
is then given by 



^5 = (1/8)|/a(0) 2 = (1/4)|/ 2 

where the physical value of / is 93 MeV. 
For (3 = 5.7, we find very consistent values, 

L 5 (f A ) = 1.66(21) x HT 3 
L 5 (\/ f A /m e ff ) = 1.63(28) x HT 3 



(10) 



L 5 (l/m e ff ) = 1.69(36) x 10 



-3 



(11) 

(12) 
(13) 



were the errors have been separately determined from a 
jackknife analysis. 

For [3 = 5.9, the results are not as consistent, reflecting 
the small residual slope found for the product, RfR m , 



L 5 (f A ) 

L 5 (yJf A /m e ff ) 
L 5 (l/m e ff ) 



= 1.51(12) x 10~ 3 (14) 
= 1.27(14) x 10~ 3 (15) 
= 1.09(15) x 10~ 3 . (16) 



Here the value determined from \jm e J' is somewhat 
smaller than the result obtain from Ja alone and is closer 
in magnitude to the L5 values found at larger beta by the 
Alpha collaboration |8| . We have discussed other aspects 
of this comparison in a previous paper |3j. 

We conclude that our analysis of the hairpin insertion 
vertex at (3 — 5.7 and (3 — 5.9 gives a striking confir- 
mation of Sharpe's analysis Q of the chiral Lagrangian 
parameters and the relation between the chiral slopes of 
f A and \/m e J f . 
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